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5!!^ '. Abstract 

^! The current status of big bang nucleosynthesis is reviewed. Particular attention is given to the 
O ' degree at which the theory is consistent with the observation of the light element abundances. 



The observational information at hand on the abundances of the four light element iso- 
topes: D, ^He, ^He, and ''Li can be characterized as being either reasonably certain (for "^He 
and ^Li); getting certain (for D); or uncertain (for ^He). In the cases of ^He and ^Li, there is 
^ ■ a growing wealth of data from extragalactic HII regions^-* on ^He and from the surfaces of old 
Population II halo stars^-* on ^Li. Indeed, we are rapidly approaching the point where our un- 
certainty in the abundances of these isotopes is dominated by systematic rather than statistical 
uncertainties. For D, though we have good solar and interstellar medium (ISM) data^'', the con- 
nection to a primordial abundance through galactic chemical evolution introduces substantial 
uncertainties. However, recent observations of D in quasar absorption systems are beginning 
to to yield a more coherent picture for what may be the true primordial abundance of D. In 
the uncertain category, I would place ^He. As for D, we have data in the ISM and there are 
also ^He abundances in planetary nebulae, however in this case, not only our we hampered by 
our lack of knowledge concerning galactic chemical evolution, but also by our uncertainties in 
the production of ^He in low mass stars. 

Overall, there is certainly a broad agreement between the predicted abundances of the light 
element isotopes from big bang nucleosynthesis^^ and the abundances inferred from observations 
which span nearly ten orders of magnitude. Indeed, the standard model theory really only 
contains one free parameter, namely, the baryon to photon ratio, rj. Agreement between theory 



and the observation of all of the four isotopes occurs when 7710 = lO^^rj is in the range, 2.8 
- 4. However, at closer level, the above range on rj is slightly high with regard to what is 
preferred by "^He, and may be too low for D and ^He when chemical evolution is combined 
with the observations. In fact, the lower bound to rj comes directly from an upper bound 
to the combined abundance of D and ^He and depends on the degree to which ^He survives 
in stellar evolution. Recent data on D and ^He may be yielding an inconsistent picture. As 
argued below, I believe that the most likely source of the problem is our treatment of ^He 
due to the considerable uncertainties in chemical and stellar evolution concerning this isotope. 
Furthermore, I will show that ^He and ^Li (abundances that we know best) are already sufficient 
in constraining the theory and that the abundance of D /H predicted from "^He and ^Li, appears 
to be consistent with the recent measurements of D/H in quasar absorption systems, though 
this requires that some of our beliefs concerning ^He need rethinking. 

Before moving to the comparison of theory and observation, it will be useful to first 
briefiy review the observational status of the four isotopes considered. More detail on each can 
be found in refs. 1-3). 

There is now a good collection of abundance information on the ^He mass fraction, Y, 
0/H, and N/H in over 50 extragalactic HII regions^"^-* . The observation of the heavy elements 
is important as the helium mass fraction observed in these HII regions has been augmented by 
some stellar processing, the degree to which is given by the oxygen and nitrogen abundances. 
In an extensive study based on the data in refs. 5) and 6), it was found®-* that the data is well 
represented by a linear correlation for Y vs. 0/H and Y vs. N/H. It is then expected that 
the primordial abundance of ^He can be determined from the intercept of that relation. The 
overall result of that analysis indicated a primordial mass fraction, Yp = 0.232 ± 0.003. In ref. 
9), the stability of this fit was verified by a Monte-Carlo analysis showing that the fits were 
not overly sensitive to any particular HII region. In addition, the data from ref. 7) was also 
included, yielding a ^He mass fraction^^ 

Fp = 0.234 ± 0.003 ± 0.005 (1) 

The second uncertainty is an estimate of the systematic uncertainty in the abundance deter- 
mination. Though the systematic uncertainty may be somewhat larger^''^ its precise value will 
be superfiuous to the discussion below. 

As I have said above, I also believe that the ^Li abundance is reasonably well known. In 
old, hot, population-II stars, ''Li is found to have a very nearly uniform abundance. For stars 
with a surface temperature T > 5500 K and a metallicity less than about l/20th solar (so 
that effects such as stellar convection may not be important), the abundances show little or no 
dispersion beyond that which is consistent with the errors of individual measurements. Indeed, 
as detailed by Spite^-* , much of the work concerning ''Li has to do with the presence or absence 
of dispersion and whether or not there is in fact some tiny slope to a [Li] = log ^Li/H + 12 vs. 
T or [Li] vs. [Fe/H] relationship. There is ''Li data from nearly 100 halo stars, from a variety of 
sources. I will use the value given in ref. 11) as the best estimate for the mean ''Li abundance 
and its statistical uncertainty in halo stars 

Li/H= (1.6 ±0.1^°:™)x 10-1° (2) 

The first error is statistical, and the second is a systematic uncertainty that covers the range 
of abundances derived by various methods. The third set of errors in Eq. (^) accounts for the 




possibility that as much as half of the primordial ^Li has been destroyed in stars, and that as 
much as 30% of the observed ^Li may have been produced in cosmic ray collisions rather than 
in the Big Bang. Observations of ^Li, Be, and B help constrain the degree to which these effects 
play a role^^\ For ^Li, the uncertainties are clearly dominated by systematic effects. 

Turning, to D/H, we have three basic types of abundance information: 1) ISM data; 
2) solar system information; and perhaps 3) a primordial abundance from quasar absorption 
systems. The best measurement for ISM D/H is^^-* 

(D/H)isM = 1.60 ± 0.09+°:?^ X 10-^ (3) 

However, it is becoming apparent that this value may not be universal (or galactic as the case 
may be) and that there may be some real dispersion of D/H in the ISM^\ The solar abundance 
of D/H is inferred from two distinct measurements of ^He. The solar wind measurements of 
■^He as well as the low temperature components of step-wise heating measurements of ^He in 
meteorites yield the presolar (D + ^He) /H ratio as D was efficiently burned to ^He in the Sun's 
pre-main-sequence phase. These measurements indicate that^^'^^-* 

= (4.1±0.6±1.4) X 10~^ (4) 
The high temperature components in meteorites are believed to yield the true solar ^He/H ratio 

0fl4,15) 

j^^^j =(1.5±0.2±0.3)xl0-^ (5) 

The difference between these two abundances reveals the presolar D/H ratio, giving, 

(D/H)© ^ (2.6 ±0.6 ±1.4) X 10"^ (6) 

Finally, there are the recent measurements of D/H in quasar absorption systems. The first of 
these measurements^^-* indicated a rather high D/H ratio, D/H ^ 1.9 - 2.5 xlO""^. However, a 
reported measurement^''^ of D/H in a second system seemed to show a very different abundance, 
D/H ^1-2 xlO~^. Most recently, a new observation^^) of the high D/H absorber was made 
resolving it into two components. The weighted average of these two components indicates 
that D/H = (1.9 ± 0.4) x 10~^ in these systems, again calling for a high primordial D/H. It is 
probably premature to use this value as the primordial D/H abundance in an analysis of big 
bang nucleosynthesis, but it is certainly encouraging that future observations may soon yield 
a firm value for D/H. It is however important to note that there does seem to be a trend that 
over the history of the Galaxy, the D/H ratio is decreasing, something we expect from galactic 
chemical evolution. Of course the total amount of deuterium astration that has occurred is still 
uncertain, and model dependent. 

There are also several types of ^He measurements. As noted above, meteoritic extractions 
yield a presolar value for ^He/H as given in Eq. (^. In addition, there are several ISM 
measurements of ^He in galactic HII regions^^-* which also show a wide dispersion 



^He^ 
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1 - 5 X 10"^ (7) 



ISM 



Finally there are observations of ^He in planetary nebulae^"-* which show a very high ^He 
abundance of ^He/H ~ 10"^. 




Each of the hght element isotopes can be made consistent with theory for a specific range 
in rj. Overall consistency of course requires that the range in 77 agree among all four light 
elements. ^He (together with D) has stood out in its importance for BBN, because it provided 
a (relatively large) lower limit for the baryon-to-photon ratio^^\ rjiQ > 2.8. This limit for a long 
time was seen to be essential because it provided the only means for bounding rj from below 
and in effect allows one to set an upper limit on the number of neutrino flavors^^-*, N^, as well 
as other constraints on particle physics properties. That is, the upper bound to Ni, is strongly 
dependent on the lower bound to rj. This is easy to see: for lower rj, the ^He abundance drops, 
allowing for a larger N^, which would raise the ^He abundance. However, for 77 < 4 x 10~^^, 
corresponding to Qh"^ ~ .001 — .002 which is not too different from galactic mass densities, 
there is no bound whatsoever on Nj^'^\ Of course, with the improved data on ''Li, we do have 
lower bounds on rj which exceed 10~^°. 

In ref. 21), it was argued that since stars (even massive stars) do not destroy ^He in its 
entirety, we can obtain a bound on 77 from an upper bound to the solar D and ^He abundances. 
One can in fact limit2i'24) the sum of primordial D and ^He by applying the expression below 
at t = 

(8) 

In d^), (yf3 is the fraction of a star's initial D and ^He which survives as ^He. For (^3 > 0.25 as 
suggested by stellar models, and using the solar data on D/H and ^He/H, one finds r^io > 2.8. 
This argument has been improved recently^^-* ultimately leading to a stronger limit^^-* r^io > 3.8 
and a best estimate r^io = 6.6 ± 1.4. The problem with this bound, is that it seems to indicate 
an inconsistency most notably in the high ^He mass fraction predicted at the large value of 
7]. It has been speculated that the cause may be underestimated systematic uncertainties in 
the ^He abundance, a problem with chemical evolution, or even a tau-neutrino mass thereby 
lowering A^,^ from 3 to 2. Indeed, at the large value of 77, if N^, is allowed to be adjusted, a value 
around 2 is needed to match the ^He mass fraction of 0.234. Even at 77, around 7710 = 3, the 
preferred value for A^j, would be well below 3^'^^-'. 

The limit 7710 > 2.8(3.8) derived using (H) is really a one shot approximation. Namely, it is 
assumed that material passes through a star no more than once. ( Although even the stochastic 
approach used in ref. 28) could only lower the bound from 3.8 to about 3.5 when assuming as 
always that > 0.25). To determine whether or not the solar (and present) values of D/H 
and ^He/H can be matched it is necessary to consider models of galactic chemical evolution^^\ 
In the absence of stellar ^He production, particularly by low mass stars, it was shown^'^^ that 
there are indeed suitable choices for a star formation rate, and an initial mass function, to: 
1) match the D/H evolution from a primordial value (D/H)p = 7.5 x 10~^, corresponding to 
rjiQ = 3, through the solar and ISM abundances; while 2) at the same time keeping the ^He/H 
evolution relatively flat so as not to overproduce ^He at the solar and present epochs. This 
was achieved for < 0.3. Even for g^ ~ 0.7, the present ^He/H could be matched, though the 
solar value was found to be a factor of 2 too high. For (D/H)p ~ 2 x 10~^, corresponding to 
7710 ~ 1.7, though models could be found which destroy D sufficiently, overproduction of ^He 
occurred unless g^ was tuned down to about 0.1. 

In the context of models of galactic chemical evolution, there is however, little justification 
a priori, for neglecting the production of '^He in low mass stars. Indeed, stellar models predict 
that considerable amounts of ^He are produced in stars between 1 and 3 M©. For M < 8M0, 
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Figure 1: The evolution of D and ^He with time. 



15 



Iben and Truran'^^^ calculate 
('He/H), 



1.8 X 10" 



Mrr 







M 



+ 0.7 



(D+ 3He)/H 



(9) 



so that at r]io = 3, and ((D + ^He) /H)j = 9 x 10"^ 5(3(1Mq) = 2.7! It should be emphasized that 
this prediction is in fact consistent with the observation of high ^He/H in planetary nebulae^"-'. 

Generally, implementation of the ^He yield in Eq. in chemical evolution models, 
leads to an overproduction of ^He/H particularly at the solar epoch^^'^^-* . In Figure 1, the 
evolution of D/H and ^He/H is shown as a function of time from refs. 14,32). The solid curves 
show the evolution in a simple model of galactic chemical evolution with a star formation rate 
proportional to the gas density and a power law IMF (see ref. 32) for details). The model was 
chosen to fit the observed deuterium abundances. However, as one can plainly see, ^He is grossly 
overproduced (the deuterium data is represented by squares and ^He by circles). Depending 
on the particular model chosen, it may be possible to come close to at least the upper end of 
the range of the '^He/H observed in galactic HII regions^^-', although, the solar value is missed 
by many standard deviations. 

The overproduction of ^He relative to the solar meteoritic value seems to be a generic 
feature of chemical evolution models when ^He production in low mass stars is included. In 
ref. 14), a more extreme model of galactic chemical evolution was tested. There, it was 
assumed that the initial mass function was time dependent in such a way so as to favor massive 
stars early on (during the first two Gyr of the galaxy). Massive stars are preferential from 
the point of view of destroying ^He. However, massive stars are also proficient producers of 
heavy elements and in order to keep the metallicity of the disk down to acceptable levels, 
supernovae driven outflow was also included. The degree of outflow was limited roughly by the 
observed metallicity in the intergalactic gas in clusters of galaxies. One further assumption was 
necessary; we allowed the massive stars to lose their ^He depleted hydrogen envelopes prior to 



explosion. Thus only the heavier elements were expulsed from the galaxy. With all of these 
(semi-defensible) assumptions, ^He was still overproduced at the solar epoch as shown by the 
dashed curve in Figure 1. Though there certainly is an improvement in the evolution of ^He, 
without reducing the yields of low mass stars, it is hard to envision much further reduction in 
the solar ^He predicted by these models. The only conclusion that we can make at this point is 
that there is clearly something wrong with our understanding of ^He in terms of either chemical 
evolution, stellar evolution or perhaps even the observational data. 

Given the magnitude of the problems concerning ^He, it would seem unwise to make any 
strong conclusion regarding big bang nucleosynthesis which is based on ^He. Perhaps as well 
some caution is deserved with regard to the recent D/H measurements, although if the present 
trend continues and is verified in several different quasar absorption systems, then D/H will 
certainly become our best measure for the baryon-to-photon ratio. Given the current situation 
however, it makes sense to take a step back and perform an analysis of big bang nucleosynthesis 
in terms of the element isotopes that are best understood, namely, ^He and ^Li. 

Monte Carlo techniques are proving to be a useful form of analysis regarding big bang 
nucleosynthesis^'^'^^'' . In ref. 36), we performed just such an analysis using only ^He and ^Li. 
It should be noted that in principle, two elements should be sufficient for constraining the one 
parameter {rj) theory of BBN. We begin by establishing likelihood functions for the theory and 
observations. For example for ^He, the theoretical likelihood function takes the form 

Lbb^Y, Fbbn) = e-(^-^BBN(.))V2.f (10) 

where YsBNiv) is the central value for the ^He mass fraction produced in the big bang as 
predicted by the theory at a given value of r], and ai is the uncertainty in that value derived 
from the Monte Carlo calculations^^-* and is a measure of the theoretical uncertainty in the big 
bang calculation. Similarly one can write down an expression for the observational likelihood 
function. In this case we have two sources of errors as discussed above, a statistical uncertainty, 
(72 and a systematic uncertainty, cXsys. Here, I will assume that the systematic error is described 
by a top hat distribution^ The convolution of the top hat distribution and the Gaussian (to 
describe the statistical errors in the observations) results in the difference of two error functions 

where in this case, Yq is the observed (or observationally determined) value for the ^He mass 
fraction. (Had I used a Gaussian to describe the systematic uncertainty, the convolution of 
two Gaussians leads to a Gaussian, and the likelihood function (0) would have taken a form 
similar to that in ([T0|) . 

A total likelihood function for each value of t]iq is derived by convolving the theoretical 
and observational distributions, which for ^He is given by 

L'^^'toUv) = J dYLBBN {Y, Ybb^ (v)) Lo{Y, Fq) (12) 

An analogous calculation is performed^^-* for ^Li. The resulting likelihood functions from the 
observed abundances given in Eqs. (|l|) and (^ is shown (unnormalized) in Figure 2. As one can 
see there is very good agreement between ^He and ^Li in vicinity of r/io — 1.8. 

The combined likelihood, for fitting both elements simultaneously, is given by the product 
of the two functions in Figure 0, and is shown in figure ^ From Figure it is clear that "^He 
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Figure 2: Likelihood distribution for each of ^He and ^Li, shown as a function of rj. The 
one-peak structure of the ^He curve corresponds to its monotonic increase with rj, while the 
two-peaks for ^Li arise from its passing through a minimum. 
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Figure 3: Combined likehhood for simultaneously fitting ^He and ^Li, as a function of 77. 
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Figure 4: D/H and ^He/H as a function of rjio from BBN along with the one a uncertainty from 
Monte Carlo calculations^^-* . Also shown are the values (demarcated by rectangles) of D/H and 
■^He/H consistent with 68% (dashed) and 95% CL (dotted) likelihood values for rjio- 

overlaps the lower (in rj) ^Li peak, and so one expects that there will be concordance, in an 
allowed range of t] given by the overlap region. This is what one finds in figure ^, which does 
show concordance, and gives an allowed (95% CL) range of 1.4 < rjio < 3.8. Note that the 
likelihood functions shown in Figures |] and ^ are not normalized to unity. An t] dependent 
normalization has however been included. Any further normalization would have no effect on 
the predicted range for rj. 

Thus, we can conclude that the abundances of "^He and ^Li are consistent, and select 
an rjiQ range which overlaps with (at the 95% CL) the longstanding favorite range around 
rjio = 3. Furthermore, by finding concordance using only ^He and ^Li, we deduce that if there 
is problem with BBN, it must arise from D and '^He and is thus tied to chemical evolution or the 
stellar evolution of ^He. The most model-independent conclusion is that standard BBN with 
Ni, = 3 is not in jeopardy, but there may be problems with our detailed understanding of D 
and particularly ^He chemical evolution. It is interesting to note that the central (and strongly) 
peaked value of rjiQ determined from the combined ^He and^Li likelihoods is at rjio = 1-8- The 
corresponding value of D/H is 1.8 x 10"'^ very close to the value of D/H in quasar absorbers^^'^^-* . 

Since D and ^He are monotonic functions of r], a prediction for t], based on ^He and ^Li, 
can be turned into a prediction for D and ^He. In Figure , the abundances of D and ^He as a 
function of rjio are shown along with the one a uncertainty in the calculations from the Monte 
Carlo results^s). The 68% (dashed) and 95% CL (dotted) ranges for D and ^He as given by our 
likelihood analysis above are shown by a set of rectangles. The corresponding 95% CL ranges 
are D/H = (5.5 - 27) x 10"^ and and ^He/H = (1.4 - 2.7) x 10"^ 

In summary, I would assert that one can only conclude that the present data on the 
abundances of the light element isotopes is consistent with the standard model of big bang 
nucleosynthesis. Using the the isotopes with the best data, '^He and ^Li, it is possible to 
constrain the theory, and obtain a best value for the baryon-to-photon ratio of rjio = 1-8 with 



a 95% CL range of 1.4 to 3.8. This is a rather low value and corresponds to a baryon density 
0.005 < Qsh^ < 0.014, and would suggest that much of the galactic dark matter is non- 
baryonic'^^-* . These predictions are in addition consistent with recent observations of D/H in 
quasar absorption systems. Difficulty remains however, in matching the solar '^He abundance, 
suggesting a problem with our current understanding of galactic chemical evolution or the 
stellar evolution of low mass stars as they pertain to ^He. 



Acknowledgments 

I would like to thank M. Casse, B. Fields, S. Scully, D. Schramm, G. Steigman, J. Truran, and 
E. Vangioni-Flam for enjoyable collaborations. This work was supported in part by DOE grant 
DE-FG02-94ER-40823. 

References 

[1] D. Kunth, these proceedings 
[2] F. Spite, these proceedings 
[3] R. Ferlet, these proceeings 

[4] T.R Walker, G. Steigman, D.N. Schramm, K.A. Olive and K. Kang, Ap. J. 376 (1991) 
51. 

[5] B.E.J. Pagel, E.A. Simonson, R.J. Terlevich and M. Edmunds, MNRAS 255 (1992) 325. 

[6] E. Skillman et al., Ap.J. Lett, (in preparation) 1995. 

[7] Y.I. Izatov, T.X. Thuan, and V.A. Lipovetsky, Ap.J. 435 (1994) 647. 

[8] K.A. Olive and G. Steigman, Ap.J. Supp. 97 (1995) 49. 

[9] K.A. Olive and S.T. Scully, IJMPA (1995) in press. 

[10] C.J. Copi, D.N. Schramm, and M.S. Turner, Science, 267 (1995) 192; D. Sasselov and D. 
Goldwirth, Ap.J. 444 (1995) L5. 

[11] R Molaro, F. Rrimas, and R Bonifacio, A.A. 295 (1995) L47. 

[12] T.R Walker, G. Steigman, D.N. Schramm, K.A. Olive and B. Fields, Ap.J. 413 (1993) 
562; K.A. Olive, and D.N. Schramm, Nature 360 (1993) 439; G. Steigman, B. Fields, K.A. 
Olive, D.N. Schramm, and T.R Walker, Ap.J. 415 (1993) L35. 

[13] J.L. Linsky, et al., Ap.J. 402 (1993) 694; J.L. Linsky, et al., Ap.J. 451 (1995) 335. 

[14] S.T. Scully, M. Casse, K.A. Olive, D.N. Schramm, J.W. Truran, and E. Vangioni-Flam, 
|astro-ph/050"8086| , Ap.J. (1995) in press. 

[15] J. Geiss, in Origin and Evolution of the Elements, eds. N. Prantzos, E. Vangioni-Flam, 
and M. Casse (Cambridge: Cambridge University Press, 1993), p. 89 

[16] R.F. Carswell, M. Ranch, R.J. Weymann, A.J. Cooke, and J.K. Webb, MNRAS 268 (1994) 
LI; A. Songaila, L.L. Cowie, C. Hogan, and M. Rugers, Nature 368 (1994) 599. 



[17] D. Tytler and X.-M. Fan, BAAS 26 (1995) 1424. 
[18] M. Rugers and C.J. Hogan, preprint 1995. 

[19] D.S. Balser, T.M. Bania, C.J. Brockway, R.T. Rood, and T.L. Wilson, Ap.J. 430 (1994) 
667. 

[20] R.T. Rood, T.M. Bania, and T.L. Wilson, Nature 355 (1992) 618; R.T. Rood, T.M. Bania, 
T.L. Wilson, and D.S. Balser, 1995, in the Light Element Abundances, Proceedings of the 
ESO/EIPC Workshop, ed. P. Crane, (Berlin: Springer), p. 201. 

[21] J. Yang, M.S. Turner, G. Steigman, D.N. Schramm, and K.A. Olive, Ap.J. 281 (1984) 
493. 

[22] G. Steigman, D.N. Schramm, and J. Gunn, Phys. Lett. B66 (1977) 202. 

[23] K.A. Olive, D.N. Schramm, G. Steigman, M.S. Turner, and J. Yang, Ap.J. 246 (1981) 
557. 

[24] D. Black, Nature Physical Sci, 24 (1971) 148. 

[25] G. Steigman and M. Tosi, Ap.J. 401 (1992) 15; G. Steigman and M. Tosi, Ap.J. 453 
(1995) 173. 

[26] N. Hata, R. J. Scherrer, G. Steigman, D. Thomas, T. P. Walker, S. Bludman and P. 
Langacker, Phys. Rev. Lett. 75 (1995) 3977. 

[27] K.A. Olive and G. Steigman, Phys. Lett. B354 (1995) 357. 

[28] C.J. Copi, D.N. Schramm, and M.S. Turner, Ap.J. 455 (1995) L95. 

[29] B.M. Tinsley, Fund Cosm Phys, 5 (1980) 287. 

[30] E. Vangioni-Flam, K.A. Olive, and N. Prantzos, Ap.J. 427 (1994) 618. 
[31] I. Iben, and J.W. Truran, Ap.J. 220 (1978) 980. 

[32] K.A. Olive, R.T. Rood, D.N. Schramm, J.W. Truran, and E. Vangioni-Flam, Ap.J. 444 
(1995) 680. 

[33] D. GaUi, F. Palla, F. Ferrini, and U. Penco, Ap.J. 443 (1995) 536; M. Tosi, G. Steigman, 
and D.S. P. Dearborn, preprint. 

[34] L.M. Krauss and P. Romanelli, ApJ, 358 (1990) 47; P.J. Kernan and L.M. Krauss, Phys. 
Rev. Lett. 72 (1994) 3309; L.M. Krauss and P.J. Kernan, Phys. Lett. B347 (1995) 347; 
M. Smith, L. Kawano, and R.A. Malaney, Ap.J. Supp., 85 (1993) 219. 



[35] N. Hata, R.J. Scherrer, G. Steigman, D. Thomas, and T.P. Walker, |astro-ph/9412087| , 
(1994). 

[36] B.D. Fields and K.A. Olive, |hep-ph/ 950834^ , Phys. Lett. B (1995) in press. 
[37] E. Vangioni-Flam and M. Casse, Ap.J. 441 (1995) 471. 



